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Abstract 

A new process associated with the nonlinear optical properties of the electromag- 
netic vacuum, as predicted by quantum electrodynamics, is described. This can be 
called photon acceleration in vacuum, and corresponds to the frequency shift that 
takes place when a given test photon interacts with an intense beam of background 
radiation. 
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1 Introduction 



It is known for a very long time that, in quantum electrodynamics (QED), 
photon-photon interactions in vacuum are possible due to the creation of vir- 
tual electron-positron pairs [1,2]. This leads to the appearance of nonlinear 
corrections of the photon dispersion relation in vacuum. In other words, the 
vacuum becomes a nonlinear medium. However, such corrections are extremely 
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weak, and for this reason the systematic study of the vacuum nonlinearity only 
started recently. In particular, proposals for the experimental observation of 
such nonlinear corrections have been published. They include photon splitting 
[3,4], second harmonic generation [6], self-focusing [5], and nonlinear wave 
mixing in microwave cavities [7]. Also, attention has been paid to collective 
photon phenomena [8], such as the electromagnetic wave collapse [9,10] as well 
as the formation of photon bullets [11] and light wedges [12]. 

In this paper, we explore another process associated with the collective pho- 
ton interactions, related to possible frequency shift of test photons immersed 
in a modulated radiation background. This new process could be called pho- 
ton acceleration in vacuum, because of its obvious analogies with the well 
documented photon acceleration processes that can occur in a plasma or in 
an optical medium [13]. Using the Heisenberg-Euler Lagrangian density, the 
appropriate dispersion relation for a photon state in vacuum in the presence 
of a generic radiation background can be found. We will consider the case 
where the background is a nearly unidirectional photon beam. The photon 
ray equations pertinent for this case is studied in detail, and the resulting fre- 
quency shift is determined. An order of magnitude estimate for the expected 
frequency shift is presented, based on analytical and numerical calculations. 



2 Dispersion relation 

The nonlinear QED effects associated with the creation of virtual electron- 
positron pairs in vacuum can be described by the Heisenberg-Euler Lagrangian 
density, which can be written as the usual classical electromagnetic Lagrangian 
density Co plus a nonlinear quantum correction 5C in the form [14] 

£ = £ + 5£ = e F + ({4F 2 + 7g 2 ), (1) 
where the quantities T and Q are determined from JF = £ /e = \(E 2 — c 2 B 2 ) 

— » — * — * — * 

and Q = c(E ■ B). Here E and B are the electric and magnetic fields, respec- 
tively. The nonlinear parameter appearing in (1) is £ = 2a 2 e 2 l h 3 /A5mfc 5 , where 
a = e 2 /2eohc ~ 1/137 is the fine structure constant. The QED corrections in 
the Lagrangian (1) are valid in the weak field, low frequency approximation, 
i.e. we require \E\ <C E crit and uj <C u e = m e c 2 /h, where E crit pa 10 16 V/cm is 
the Schwinger critical field, u e is the Compton frequency, and uj is the typical 
frequency with which the electromagnetic field changes. These approximations 
ensures that there is no appreciable pair creation due to multi-photon effects 
(as it will be exponentially suppressed for low field strengths) and that there 
are no single photons able to generate pairs from the vacuum. However, it is 
worth pointing out that we do not require the fields to be constant in time or 
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space, only slowly varying with respect to the Compton frequency [2]. As this 
is the case for almost all relevant fields, the applicability of the Lagrangian (1) 
is guaranteed for a wide variety of field configurations. The resulting Maxwell's 

equations in vacuum take the usual form, if we use the following definitions 

—»—»—»—» —»—» — » 

D = e E + P and B = jj, (H + M), where the polarization P and magneti- 
zation M result from the nonlinear QED corrections included in (1), and are 
given by 



P = 2C(4^E + AcGB) , M = -1<?C,{\TB + 7QEjc). (2) 

Starting from Maxwell's equations in vacuum, we can then establish the equa- 
tion of propagation for the electric field E in the form 



(V 2 - c~ 2 dt) E = /2 [d 2 P + c 2 V(V • P) + d t iy x M)\ , (3) 

and a similar equation for the magnetic field. If we consider a given photon 
state, with the frequency u and the wavevector k, we can derive from these 
wave equations the following nonlinear dispersion relation [3,7] 



uj = A; C (l-lA|Q| 2 ), (4) 
where Q represents the background electromagnetic field. We have \Q\ 2 = 



e 



and A = A±, with A + = 
arization states. We have 



E 2 + C 2 B 2 _ . E y _ c 2(- . B y _ 2cH . . _g) 

14^ and A_ = 8( for the two independent photon poi 
also used the unit propagation vector n = k/k. We note that Eq. (4), to first 
order in the QED corrections, can be rearranged to give the wavenumber as 
a function of frequency (both forms being equivalent). However, the current 
form will be used in Sec. Ill (see Eq. (9) and the discussion following it) as it 
is the canonical choice using the geometric optics approximation. 

It may also be argued that derivative corrections to the weak field Heisenberg- 
Euler Lagrangian (1) could be of importance in certain applications, especially 
for high-frequency field configurations where the local field approximation 
may be questionable. However, the derivative correction to the weak field 
Lagrangian (1) was given in Ref. [15], taking the form 



2e c 2 a 

oL D 



(d a F ab )(d c F%)-F ab d 2 F ab , (5) 



where F ab is the field strength tensor and we have used the summation conven- 
tion over the four-indices a, b, c, and this correction gives rise to terms in the 
dispersion relation proportional to the parameter a 2 (uu /oj e ) 2 (\E\/ E CT - lt ) 2 [16]. 
Thus, for frequencies much smaller than the Compton frequency, appropriate 
for the field configurations considered below, these derivative correction are 
smaller by a factor a(uo /uo e ) 2 1 than the QED corrections in the Lagrangian 
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(1), and they may therefore at this stage be safely neglected in comparison 
to the correction in (1) (that are proportional to a(\E\/E crit ) 2 ). However, on 
longer time-scales such corrections could in principle yield interesting disper- 
sive effects on the interaction among photons [11,12]. For a thorough discussion 
of the Euler-Heisenberg Lagrangian and its validity limits, see e.g. Ref. [17] 

Let us now consider the simplest case where the background radiation field is 
associated with a single photon state with the frequency uj' and the wavevector 
k! . Here the dispersion relation (4) takes the form 



kc 



l-h \f(k,k')\E(k 



f \ 1 2 



(6) 



— * — » — * — * 

where E(k') is the electric field of the background radiation, and f(k, k') is a 
geometric factor determined by 



f(k, k') = {2 - (n ■ e*) 2 - [n ■ (n' x e*)] 2 - 2(ra • ft')]} (7) 

— * 

involving the unit vector n defined above, as well as n' = k'/k' and e 1 = 

— » — » — » — * — * 

E(k')/E(k'). It is obvious that, for n = ft', we have f(k, k') = 0, which means 
that there is no coupling between photons belonging to the same photon state. 
In other words, the nonlinear self-interaction of a photon with itself is strictly 
forbidden (to all orders in the amplitude, see Ref. [2]). Only different states of 
the radiation background will influence the dispersion of a photon in vacuum, 
through the virtual electron-positron pair fluctuations. A test photon moving 
in a background field will generically have a non-vanishing geometric factor 
(7). In particular, a background field consisting of a plane wave will result in a 
non- vanishing geometric factor unless the test photons are propagating parallel 
to the wave vector of the background field. We emphasize that the finite value 
of the geometric factor (10) is a direct consequence of non- vanishing values of 
the relativistic invariants F and G when considering the total electromagnetic 
field, i.e. incorporating both the background fields and the test photon field 
in the calculation scheme. 

We can easily generalize the dispersion relation (6) for the case of an arbi- 
trary broadband radiation spectrum, which can be slowly modulated in space 
and time over length and time scales much longer than the periods and the 
wavelengths of the considered radiation states. The resulting expression is 



uj(f, k, t) = kc 

where we have introduced the photon occupation number, or number of pho- 



f(k,k')k'N(f,k',t)-^y 3 



(8) 
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tons per field mode, as determined by 



Ahk 

and a new coupling factor determined by A* = AHX/c. 



3 Photon dynamics 



In the geometric optics approximation, the dynamics of photons can be de- 
scribed by the ray equations, which can be stated in the following canonical 
form [18] 



df duo dk duo 

dt = ~dJ ' dt = ~W' ( } 

where the Hamiltonian uo is determined by Eq. (8). From (10) follows that 
duo/dt = duo/dt, and in accordance with the general principles of geometrical 
optics, the ray equations imply that a time-dependent medium is needed for 
photon acceleration (frequency conversion) to occur, i.e. for a pure space- 
dependence of the Q-factor in (4) the photon frequency is conserved. 

Let us consider the important particular case where the background radiation 
is dominated by a single photon beam propagating in a direction Oz. This 
means that we can make k' = k'e z . If the probe photon described by these 
equations of motion propagates at a given angle 9 with respect to the z-axis, 
we have (n ■ n') = cos 9 and (ft • e f ) = sin#cos-?/>, where ip is the second 
angle necessary to define the direction of n with respect to the background 
electric field. From (7), we obtain the geometric factor f(k,k') = f(9) = 
[2(1 — cos#) — sin 2 9]. The photon dispersion relation then becomes 



co(f, k, t) =kc[l- \\*f{9)I{r, t)] , (11) 

where the angle 9 can vary in space and time, according to the photon dy- 
namics described by Eq. (10), and the background field intensity I(f;t) is 
determined by the integral 



r -> dk' 

I(f,t) = j k'N(r,k>,t) — . (12) 

Let us assume that the beam is modulated in intensity along the propagation 
direction, and that we can write I(z — ut,f±). For a very large beam waist, 
we can neglect the dependence over the transverse direction and assume that 
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u = c. From the equations of motion (10) we can then conclude that the 
transverse photon dynamics is an invariant 



— — = —t^- = , A;_i_ = k sin 9 = constant. (13) 
at or± 

Even if this is not exact, a small variation in k± would only result in a slight 
diffraction effect, with no significant contribution to the frequency shift that 
will be examined below. In the general case, however, both k and 9 will evolve 
in space and time, as determined by the parallel equations of motion 



dz 



duo 
dk z 



k z c 
k 



l-^f(9)I(z-ct) 



and 



dk z 



(14) 
(15) 



The parallel photon motion can be more appropriately described by using a 
new variable r) — z — ct. This can be formally done by making a canonical 
transformation from (z,k z ) into a new pair of variables (r),p), through the 
generating function F(z,p,t) — (z — ct)p. Hence, we have 



OF OF 
k z = — =p , r] = — = z - ct. 
oz op 

The resulting new Hamiltonian will then be defined by 



(16) 



u'(rf,p) = u;(rj : p) + 



OF 



or, in explicit form, by 

/{i],p) = c\J k\ + p 2 



1 - yiWM 



cp, 



(17) 



(18) 



where kj_ is a constant, as stated by equation (13). The canonical equations 
for the parallel photon motion become 



and 



dr) duj' 
dt dp 

dp 



-c + 



pc 



1 - 



(19) 
(20) 
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In order to understand the implications of these equations, we will assume that 
the beam intensity is modulated over a distance A, as described explicitly by 
the function 



I(rj) = I g( V ) = | [1 ± tanhfa/A)] , (21) 

where the signs + and — can describe the rear and front of the beam , respec- 
tively. 

Before describing any explicit solution, it is however important to notice that 
the new Hamiltonian u/ is an invariant, which allows us to establish a relation 
between the values of the photon frequency at two distinct positions of its 
trajectory. Let us then consider two different positions along the same photon 
trajectory, viz. r\ = rji and r\ = rj 2 - Because of the invariance of u/, we can 
write us' = k\c{l — Sni) — cpi = k 2 c(l — 5n 2 ) — cp 2 , where we have used 
ki = k(r)i) and pi = p(r)i), for i — 0,1, and defined the nonlinear refractive 
index perturbations 5rii = |A* f(0i)I(r)i). 

Notice that, far ahead of (or far back of) the background beam we have Sn xs 
0, and deep inside the beam we attain the maximum perturbation Sn max ~ 
(A* /2)f(9)I . This expression can be explicitly written in terms of the photon 
frequencies noticing that they are related to the wavenumbers ki through 
the dispersion relation (11), which can be rewritten as 



This leads to the following relation between the initial value of the frequency 
c^o = oj(t = 0), and a subsequent value uj\ = u>(t) 



1 - Sn(t) (1 - cos - 5n ) 
^ ~ ^ (1^) l- cos 9(t)-5n(ty (23) 

This expression will be very useful to calculate the frequency shift of a test 
photon due to its interaction with the radiation background. Such a frequency 
shift can be associated with a process of photon acceleration. Notice that the 
group velocity of the photon as determined by the nonlinear dispersion rela- 
tion, and given by [3] v = nc[l — 5n(t)], is independent of the photon frequency, 
but changes with time due to the change in the angle between the photon and 
the background beam. In order to determine how this process evolves with 
time, we have to go back to the parallel photon dynamical equations. It is 
useful to rewrite them using a new time variable, defined by r = ct A*/ /2, 
which describes the photon evolution on a very slow time scale. Introducing 
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Fig. 1. (a) Photon acceleration process at the rear of an intense beam in vacuum: 
Evolution of the relative parallel wavenumber £(t)/£(0), and the relative frequency 
shift Y(t) (in bold), as a function of r/A, for no = 1, for an initial value £(0) = 5. (b) 
Photon deceleration process at the front of an intense beam in vacuum: Evolution 
of the relative parallel wavenumber £(r) /£(()), and the relative frequency shift Y(t) 
(in bold), as a function of r/A, for uq = 1, for an initial value £(0) = 5. 

the dimensionless variable £ = p/k±, we have then 



dr) 
dr 



2 i \^ + X(l-Xf 9 (r ] ) 



and 



§ = (i-*) 2 v / l + CV(^ 



(24) 



(25) 



where X = cos 6 = £/ \/\ + £ 2 . For a plausible situation, even for very intense 
background radiation intensity, such that A* Jo <^ (1 — X), we can use the 
approximate solution 



„(r) « - 2[1 : x(0)] 7- 



"MOT, 



and, using Eq. (21), we can write 



(26) 



§ « ±[1 - *(C)] 2 ^\/lTe sech 2 (—u T/ A). 



(27) 



This approximate description of the photon acceleration process occurring in 
vacuum, on a very slow time-scale defined by r, is completed by Eq. (23). It 
is nearly equivalent to 



Y(t) 



u{r) ^ 1-X(0) 
w(0) ~ 1 -X(t)' 



(28) 
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The results are illustrated in Fig. la , where the values of £(r) and Y(t) are 
represented for photons with initial values ((0) = 5. Here the photon is inter- 
acting with the rear of the background beam, and we can see that a noticeable 
up-shift of the photon frequency can take place. This is an interesting result, 
given the extremely small value of the nonlinear QED term responsible for 
the interaction. The physical meaning is, however, very clear: the very weak 
nonlinear force acting on the test photon and due to the gradient of the beam 
intensity (represented by the derivative g'(i])) acts on the photon over ex- 
tremely large distances, because it travels with a parallel velocity nearly equal 
to the light speed c, thus leading to a non- negligible frequency up-shift. No- 
tice that the approximate solution (27) gives an overestimate of the relative 
velocity uq, and therefore underestimates the up-shift. The frequency shifts 
tend to increase for lower values of C(0)- 

The effect of a frequency down-shift can occur at the front of the beam, as 
illustrated in Fig. lb. We see that, in this case, the photon frequency deceases, 
and can eventually tend to zero, for lower initial values of (. 



4 Conclusions 

We have described here a new process associated with the nonlinear optical 
properties of the electromagnetic vacuum, as predicted by quantum electrody- 
namics. This can be called photon acceleration in vacuum, and can take place 
when a given test photon interacts with a modulated background radiation. 
This process is similar to that already observed in laboratory experiments, 
using relativistic ionization fronts [19]. Here, however, there is no material 
supporting the interaction, except for the QED nonlinearities of vacuum. The 
ionization front is replaced by a trail of virtual electron-positron pairs created 
by the intense background radiation in vacuum. In the present case we have a 
kind of virtual ionization front, which is responsible for the photon frequency 
shift. 

The case when the background is made of a single intense beam, propagating 
along a given direction, was considered in detail. It was shown that, over 
long interaction distances, the direction of propagation and the frequency of 
test photon can be significantly changed. A frequency up-shift occurs at the 
rear of the beam, and a down-shift can take place at the front. A similar 
process can occur when a photon interacts with an acoustic perturbation of 
the background, and will be the object of a forthcoming publication. In this 
case, a stronger effect is to be expected, because the geometric factor f(6) will 
be larger than in the present case of nearly parallel photon-photon interaction. 

In the present work only individual photon dynamics was considered. But the 
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test pulse can belong to the radiation background itself, which means that the 
process described here will lead to spectral changes of the radiation spectrum 
in vacuum. This spectral evolution can be conveniently treated with the help 
of photon kinetics [13], to be analyzed in the future. 
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